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ABSTRACT 

We simulate clouds in the Galactic Centre (GC) crossing over the black hole in parts 
and present this as a possible formation mechanism for the observed stellar disks in 
the GC through the redistribution of angular momentum by colliding material with 
opposite angular momentum. A parameter study using six high resolution simulations 
of an isothermal cloud of constant density falling onto the black hole and crossing over 
it in parts demonstrates that this mechanism is able to reproduce the observed disk 
properties in the GC. The evolution of the ensuing accretion disks is highly non- linear 
with the redistribution of the angular momentum through dissipative processes being 
a dominant effect. We analyse the resulting Toomre unstable, eccentric gaseous disk 
and show that this already yields a good comparison with the observed stellar disk 
size and eccentricity in the GC. The best simulation results in an outer radius of 1 pc, 
a mass of 10^ Mq and an eccentricity of 0.24 for the Toomre unstable disk, which 
compares well with the observations. 
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1 INTRODUCTION 



Observations reveal an interesting feature in our Galactic 
Centre (GC): two sub-p arsec-scale rings of younR stars near 
the radi o source SgrA* (jCenzel et al"]|2003l : IPaumard et al.l 
1.2006; B artko et al.ll2009r ) which cannot be explained by nor- 
mal means of star formation due to the hostility of the 
environment. Tidal forces would disrupt typical molecular 
clouds in the vicinity of the central black hole, preventing 
their condensation into stars. There are around 100 young 
and massive stars distributed in a warped clockwise rotating 
disk and a second inclined counter-clockwise rotating disk. 
The outer edge of the system is at around 0.5 pc and the 
inner edge at around 0.04 pc. The mean eccentricity of the 
clockwise rotating system is measured to be 0.36 ± 0.06. 

An alternative and more plausible formation scenario 
postulates that the stellar rings formed by fragmentation of 
self-gravitating, eccen tric accretion disks (|Rice et al ] |2005l : 
[Alexander et al.ll2008l and references therein) . Previous stud- 
ies concentrated mainly on the fragmentation of an al- 
ready existing accretion disk. The question how these disks 
formed i n the first place has only recently be en discussed , 
e.g. by (IWardle fc Yuscf-Zadch 2008; Mapcl h et al] l2008l : 
iHobbs fc Navakshin. ,200a : .Bonnell fc Rica ,20di '). In their 
models the accretion disks are built up by the rapid deposit 
of gas around the central black hole through the infall and 



tidal disruption of a gas cloud. However, here the problem 
arises that a gas cloud that most likely formed far from the 
GC where tidal forces are inefficient must be placed on an 
orbit with a close passage around the black hole. 

Furthermore, the cloud should not be disrupted by 
internal st ar formation proces s es before encountering the 
blac k hole. iMapelh et al] ||2008| ) , iHobbs fc NavakshinI (|2009D 
and iBonnell fc Ricd ~ 20081 ) studied the formation of a 
disk by a cloud or clouds passing the central black hole 
at some distance large r than the cloud radius. However, 
IWardle fc Yusef-Zadehl (|2008l ) showed that it is very un- 
likely for a cloud to be captured by the black hole with- 
out engulfing it during the encounter, as only a small set 
of initial parameters would lead to such an event. Thus 
IWardle fc Yusef-Zadehl (|2008l ) proposed a model in which 
the infalling cloud covers the black hole in parts during its 
passage. This process efficiently redistributes angular mo- 
mentum and dissipates kinetic energy through the collision 
downstream of material with opposite angular momentum 
streaming around the black hole from both sides, finally re- 
sulting in a compact accretion disk. 

In this paper we study the direct encounter of a gas 
cloud with the central black hole in detail using high- 
resolution numerical simulations. We show that gas disks 
similar to the observed stellar GC disk can be formed when 
the infalling clouds cover the black hole in parts. This paper 
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is structured as follows. In Section[2]we present a motivation 
for our model. In Section |3] the numerical method and the 
initial setup of the model are presented. The main results 
as a function of varying impact parameters and initial cloud 
velocities are presented in Section |4] Finally, we summarize 
and discuss our findings in Section [5] 



2 MOTIVATION 

In order to motivate our model we start by considering a 
cloud on an arbitrary Keplerian orbit around the black hole, 
where the cloud settles into a circular orbit by dissipating 
kinetic energy into thermal energy. The largest cloud which 
could form a disk similar to the observed stellar disk of ra- 
dius rdisk ~ 0.5 pc, which at the same time does not cross 
the black hole during the passage, would be a cloud of ex- 
actly the stellar disk radius of 0.5 pc on an initially circular 
orbit. All clouds with a larger radius will cross over the black 
hole in parts during the passage. The minimum distance to 
the black hole must be smaller than Vdisk for all non-circular 
orbits if we want to build up a disk of size Vdisk by energy dis- 
sipation because of angular momentum conservation. Thus 
the circular case really gives the upper limit for the size of a 
cloud that does not touch the black hole during the passage. 

The required mass in the gaseous disk that 
formed the observed stellar disk is around 1 0^~^ Mq 
l|Navakshin fc Cuadral l2005l : iNavakshin et al.l |2006|) . A 
cloud of radius 0.5 pc with a mass of 10* Mq already has a 
density above 10^ cm~^ which is the obse rved upper density 
limit for molecular clouds in the GC l|Giisten fc Philipa 
|2004 ). 

Furthermore, if we can only form a disk of given radius 
by a cloud losing velocity due to dissipation of energy then 
all incoming clouds must have an orbit where the shortest 
distance to the black hole is smaller than the resulting disk 
radius. In the case of angular momentum redistribution this 
is no longer necessary. Even clouds on orbits with a point 
of closest approach much larger than the wanted disk ra- 
dius can potentially form such a disk if their radius is larger 
than the minimum distance to the black hole so that angular 
momentum can be redistributed. 

Thus, we conclude that the redistribution of angular 
momentum must always be included as a mechanism for ini- 
tiating the formation of such a disk. This will also increase 
energy dissipation by the formation of colliding flows. The 
aim of this paper is to show that this mechanism can in- 
deed form a disk similar to the one observed in the GC, 
wh ich has alread y been sho wn for pure ene r gy di ssipation 
bv lMapeUi et al.l f2008i l and lBonnell Sz Ricd l|2008l ). 



3 MODEL AND NUMERICAL METHOD 

The hydrodynamical evolution of the impact of the gas cloud 
with the central black hole is studied using simulation runs 
performed with the N-body Smoothed Part icle Hydrody- 
namics (SPH) Code Gadget2 l|Springelll2005l ). The simula- 
tions are run with a total number of 10® SPH particles and a 
softening length of e = 10~^ pc. The number of neighbours 
is set to rincigh ~ 50 ± 5. The black hole is included as a 
static potential of a point mass of Mbh = 3.5 X 10® M0 
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ID 7j t)c [km/s] b [pc] jspociflc [jf^] 



101 


0.28 


0.29 


50 


1 


50 


102 


0.57 


0.41 


50 


2 


100 


vol 


0.85 


0.3 


30 


3 


90 


V02 


0.85 


0.5 


50 


3 


150 


V03 


0.85 


0.8 


80 


3 


240 


GDI 


0.57 


0.98 


120 


2 


240 



Table 1. Summary of initial conditions for our simulations. The 
cloud velocity is varied for simulations VOl, V02 and V03 and 
the impact parameter for simulations 101, 102 and V02. COl is 
a comparison simulation with the same initial specific angular 
momentum (jspcciflc) V03. 

("Cen zel et al]l2003h placed at the origin of the coordinate 
system. 

We adopt an isothermal equation of state with a typi- 
cal c loud temperature of Tdoud = 50 K ()Giisten fc Philipg 
|2004 ). This simplification should not strongly affect our re- 
sults, i.e. the formation of the accretion disk, as the orbital 
velocities (uorbit '~ 100 km/s) are in general much larger 
than the cloud sound speed {cs ~ 0.45 km/s). Heating and 
cooling processes in the disk will however play an impor- 
tant role when investig ating its fragmentation and condensa- 
tion into stars (see e.g. iBonnell fc Ricell2008l ). We postpone 
this question to the next paper in this series where we will 
study in detail the thermodynamics of the accretion disk, 
viscous heating, radiative cooling, star formation and as- 
sociat ed feedback processes, as detailed in Johansso n et al.l 
(l2009li . 

In this first study we start with a spherical and homo- 

J eneous cloud of radius 3 .5 pc which is typical for the GC 
Mivazaki fc Tsuboil [ioool . derived from CS(l-O) radio sur- 
veys) . The H2 gas density is 10* cm~^ (|Giisten fc PhilippI 
|2004| . taken from high-resolution surveys of CS and CO in 
the GC), leading to a cloud mass of 8.81 x 10* Mq. The 
SPH particle mass is then mspH = 8.81 x 10~^ Mq and 
the corresponding minimum mass that can be resolved is 
mmin = nncigh X mspH = 4.4 Mq. 

In order for the simulations not to become too time- 
consuming as a result of very small particle time steps in 
the vicinity of the black hole we define an accretion radius 
^acc within which all SPH particles are considered to be 
accreted by the black hole and are removed. We want this 
accretion radius to be larger than the minimum smoothing 
length of 10~^ pc and smaller than the observed inner disk 
edge of 4 X 10~^ pc. Test simulations showed that a value 
of race = 2 X 10~^ pc captures the most relevant accretion 
physics and disk structure while at the same time allowing us 
to simulate the encounter until the whole cloud has passed 
the black hole and an accretion disk has formed. 

We do not include black hole growth in order to speed 
up the simulations. This growth is anyhow negligible since 
at most only around one percent of the black hole mass 
was accreted at the end. Energetic feedback from the ac- 
creting black hole might however be important and will 
be investigated in a subsequent paper. The gravitational 
forces are purely Newtonian, relativistic effects (frame drag- 
ging) would only become important at radii smaller than 
r = 10~* pc for a Af = 3.5 x 10® Mq black hole which is 
beyond our accretion radius. The Gadget Code treats grav- 
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Figure 1. Logarithmic density in in the xy-plane with z=0 pc at different times for simulation 102. Velocity is given in km/s with 
the unit length of the arrow corresponding to 148 km/s. Material with opposite angular momentum collides, a density wave builds up 
and spirals into the forming disk. 



itational forces as Newtonian up to 2.8 times tlie smootliing 
lengtli, wliicli in any case is also beyond our accretion radius 
for the minimum smoothing length. 

Initially, the centre of mass of the cloud is placed at 
a distance of 5 pc from the origin/black hole on the x-axis 
(the direction of motion) and a offset 6 on the y-axis. We 



have chosen this relatively small distance in order to prevent 
the cloud from collapsing before it reaches the black hole. 
The free fall time of a cloud with n = 10* cm~"^ is only 
around 0.36 Myrs. For an average infall velocity of 50 km/s 
the cloud's centre of mass should reach the black hole after 
roughly 0.1 Myrs. 
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Figure 2. A zoom-out of the logarithmic density plot from[lJ' in 
in the xy-plane with z=0 pc. Velocity is given in km/s with 
the unit length of the arrow corresponding to 199 km/s. 



The evolution is characterized by the fraction where 
R is the cloud radius and b is the initial cloud's centre of 
mass offset on the y-axis. is always smaller than one in 
our case of a cloud crossing over the black hole. 

The initial cloud velocity Vc in x-direction is a typ- 
ical value, taken from the observations presented by 
iMivazaki fc Tsuboil (|2000l '). The evolution is then character- 
ized by the ratio ^ where Vc is the initial cloud velocity 
and vl — '^^'^^>'f>- is the black hole's escape velocity at dis- 
tance 6. Table [1] summarizes the initial conditions we used 
for our parameter study. Three simulations were computed 
with fixed initial cloud velocities of 50 km/s and varying 
offsets 6, set to 1, 2 and 3 pc, respectively. In addition we 
did two simulations with b fixed to 3 pc and varying initial 
cloud velocities, set to 30 and 80 km/s, respectively. The 
comparison simulation COl has the same initial specific an- 
gular momentum as simulation V03 but a smaller impact 
parameter and higher velocity. 

The simulations were typically run for an evolutionary 
time of 0.25 Myrs which corresponds roughly to the situation 
when all parts of the cloud have crossed the black hole. Only 
simulation 101 was stopped earlier, that is after 0.2 Myrs. 
At that time, due to the small impact parameter, the time 
step became very small making the simulation very time- 
consuming and too expensive to continue. In this paper we 
are only interested in the resulting gaseous disk properties 
and defer the study of star formation to a follow-up paper. 



ID 


A/acc,exp [10* Mq] 


A/acc,sim [10* Mq] 


101 


0.1 


4.33 


102 


0.08 


1.22 


vol 


0.13 


0.68 


V02 


0.05 


0.26 


V03 


0.03 


0.16 


COl 


0.03 


1.12 



Table 2. The total mass within the accretion radius iWacc ,exp as 
expected theoretically, assuming angular momentum conservation 
compared with the result A/acc.sim from the numerical simula- 
tions. 

4 RESULTS 

4.1 Density Evolution 

Figure [1] shows the time evolution of the cloud for our stan- 
dard model 102. Density maps and the velocity field in the 
xy-midplane with z=0 are plotted. After starting the simu- 
lation, the cloud approaches the black hole from the right 
side along the x-axis. Due to tidal forces the parts closest to 
the black hole quickly start to form a finger-like extension 
stretching towards the black hole located at the coordinate 
origin as seen in Figure [T^. 

This material passes the black hole on orbits corre- 
sponding to its initial angular momentum. 8000 yrs later 
(Figure [TJs), this material collides with gas from the oppo- 
site side. A density wave forms (green area) which thickens 
and moves around the disk as seen in the figures [T]:,d,e. 
The simulation was stopped after 0.25 Myrs when the cloud 
had passed around the black hole. An inner, relaxed equi- 
librium disk has formed at that time while the outer region 
represents a 1-armed infalling spiral in the disk plane. This 
evolutionary state is shown in Figure [Tf and Figure [5] Note 
also the narrow yellow high-density line along the negative 
X-axis in Figure [TJ; and Figure [TJl which results from the col- 
limation of material, initially located above and below the 
z=0 plane which after passing the black hole is redirected 
towards the negative x-axis. 

In all simulations, the SPH particles in the inner 0.1 pc 
reach the minimum smoothing length so that we are limited 
by numerics in this region and the physics is not resolved 
anymore due to strong smoothing over the forces and thus 
we cannot make any reliable predictions for the size of the 
inner disk radius. 

4.2 Disk Surface Density Distribution 

The resulting disk surface density distribution at the end 
of the simulations (at t = 0.25 Myrs for all simulations, 
except 101 which stopped at t = 0.2 Myrs) are shown in 
Figure [3] The surface density was calculated by integrating 
the density in the z-direction from — 10~^ to lO"'^ pc in order 
to capture only gas within the thin disk. The vertical black 
line in Figure |3] indicates the radius of 2 x 10~^ pc within 
which SPH particles were removed from the simulation and 
are considered "accreted" onto the black hole. 

The green (dashed) line shows the expected surface den- 
sity if we just redistribute material according to its initial 
angular momentum in a circular Keplerian disk. Given in 
red (solid line) is the resulting surface density from the last 
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Figure 3. Logarithmic surface density in 



'^0 



for the final time step of all simulations in red (solid line) compared with the expected 



surface density from redistributing particles in a circular Keplerian disk according to their initial angular momentum in green (dashed 
line) . 



time step of our simulations. Without angular momentum 
redistribution we would expect that the total accreted mass 
(accretion will be discussed in detail in subsection I4.3|l is 
equal to the mass within the accretion radius, given by the 
surface density distribution according to the green curve. 



In Table[2]we compare the accreted mass expected from 
this simple model to the actual accreted mass at the end of 
our simulations, which take angular momentum redistribu- 
tion realistically into account. The simulations result in val- 
ues that are a factor of 5-40 larger than naively expected 
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Figure 4. Plot of all grid cells with a Toomre parameter smaller than one. Each grid cell has a size of 0.02x0.02 



from redistributing material in a circular Keplerian disk 
according to the initial angular momentum which demon- 
strates the highly non-linear evolution of the accretion disk, 
the importance of dissipative processes and the redistribu- 
tion of angular momentum. 

To analyse disk stability we calculate the Toomre pa- 
rameter Q = on a 250x250 grid applied to an area of 
5x5 pc in the xy-plane. Figure |4] shows the resulting disks 
defined by all grid points with a Toomre parameter smaller 
than one (Q < 1). We compare the properties of the unstable 
disks defined in this way to the observed stellar disk. 

4.3 Accretion 

Accretion rates are calculated from material which falls be- 
low our accretion radius of race = 2 x 10~^ pc. This does not 
necessarily mean black hole accretion since the formation of 
a small and hot black hole accretion disk evolving viscously 
is beyond the resolution limit and outside the scope of our 
current simulation. Nevertheless, it is interesting to compare 
our accretion rates for all simulations to the Eddington rate 
of the black hole which is shown in Figure [5] The Eddington 
mass accretion rate is defined as M^dd ~ ^I^MmjUHs. with 
rUp being the proton mass, e the accretion efficiency and gt 
the Thompson scattering cross-section. From this we get an 

Eddington rate of M^dd = 0.0775 using an accretion 
efliciency of 10% and a 3.5 x lO** Mq black hole. 

Note that the y-axis is scaled differently for different 
simulations due to the large range of accretion rates we get. 
At the end of all simulations we reach very low accretion 
rates indicating that we have reached a nearly stationary 
state. The low impact parameter simulation 101 accreted al- 
most the same mass as a test case we did with -^=0. For the 
■^=0 case all cloud material has zero net angular momen- 
tum and the complete cloud is accreted in a small stream 
onto the black hole. The evolution of this test-case can be 
seen in Figure [6] together with a close-up of the accretion 
back-flow. The accretion rates of the large impact parame- 
ter simulations V02 and V03 stay below the Eddington rate 
at all times. 

A comparison of the simulations 101, 102 and V02 in 
Table [2] shows that the impact parameter plays the most 
important role in determining the final accreted mass while 
the cloud velocity has a much smaller effect (compare VOl, 
V02 and V03). In all simulations we still have accretion at 
the end of the runs, however at very low rates, thus the 
values shown in Table [5] would only increase slightly if the 
simulations were continued longer. 

4.4 Disk Properties 

Finally, we study the masses and eccentricities of our result- 
ing disks from section 14.21 Simple elliptical fits were made 
to the disks and the masses were calculated by counting gas 
which is unstable according to the Toomre criterion. The 
results are summarized in Table O 

When comparing the disk masses in the simulations 101, 
102 and V02 (fixed velocity, variable impact parameter) with 
the simulations VOl, V02 and V03 (fixed impact parameter, 
variable velocity) we note that in contrast to the accretion 
rates, for which the impact parameter was the dominating 
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Simulation A/disk [10'' Mq\ 


e 


long, short-Axis [pc] 


101 


3.54 


0.51 


1.05, 0.9 


102 


4.4 


0.41 


1.65, 1.5 


vol 


6.95 


0.52 


1.40, 1.2 


V02 


3.36 


0.33 


1.70, 1.0 


V03 


0.79 


0.24 


1.00, 0.97 


COl 


0.3 


0.54 


0.47, 0.4 


Observation 


~1.0 - 10.0 


0.34 ± 0.06 


0.5 



Table 3. Disk properties, including the mass, eccentricity and 
long, short-Axis of the disks from the last simulation time step 
(0.25 Myrs). Note that we compare properties of our gaseous disks 
to the observed stellar disk. The mass in the observations is what 
we would e xpect for a gaseous disk that formed the observ ed 
stellar disk l|Navakshin fc Cuadrall2005l : iNavakshin et al.ll2006l '). 

factor, the cloud velocity is the most important factor in 
determining the disk masses. 

Comparing V03 with COl we see that the initial specific 
angular momentum does not det ermine the final disk p rop- 
erties. According to the data from lPaumard et al.l (|2006l ) the 
stars in the observed disks have specific angular momentum 

2 

reaching from roughly 30 to 70 101 is within this range 

and vol close to it, but these simulations do not reproduce 
the properties of the observed disks as well as V03 which 
has a rather large initial specific angular momentum, which 
again demonstrates the highly nonlinear evolution of a cloud 
crossing over the black hole. 

On average we find disk eccentricities of order e = 0.4 
and an outer radius of r 1.35 pc. The observed eccen- 
tricities of the stell ar rings in the GC are e = 0.34 ± 0.06 
(|Bartko et al.ll2009l ). the estimated mass of the gaseous disk 
that f ormed the stars is expected t o be around M ~ 10**" ^ 
M0 (iNavakshin fc Cuadral l2005l : iNavakshin et all 120061 '). 
The observed outer radius of the disk is roughly r = 0.5 pc 
(|Genzel et al.|[2003 ) , with the observed disk size still increas- 
ing due to the addition of newly observed stars belonging to 
the disk structure. 

These properties are very similar to our gaseous disks, 
which have a bit too large radii but already show the same 
low eccentricity as the stellar disks and the required mass in 
gas to form the observed stellar disk. 



5 SUMMARY AND DISCUSSION 

We have performed simulations of a gas cloud colliding with 
a massive black hole with parts of the cloud crossing directly 
over the black hole during the process. Despite the fact that 
the impact velocity of the cloud and its initial distance is 
large, sub parsec-scale gas disks form due to angular mo- 
mentum redistribution and through the efficient dissipation 
of kinetic energy by the collision of gas fiowing around op- 
posite sides of the black hole. 

For non-zero impact parameters the net angular mo- 
mentum generates an eccentric gas disk with global prop- 
erties that are similar to the observed stellar rings de- 
tected in the GC. We get low eccentricities in the range 
of e = 0.24 — 0.51, disk sizes from r = 1 — 1.7 pc and masses 
of the unstable disk region from M = 0.7 - 6.95 x lO" Mq. 
The impact parameter is the dominant factor in determin- 
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Figure 5. The red (solid) line shows the accretion rate in (left y-axis). The dashed green line shows the Eddington limit of a 

3.5 X lO'' M0 black hole at 10% accretion efficiency. Also plotted is the total accreted mass in blue (dotted line) in Mq (right y-axis). 
The low impact parameter simulations accrete at super-Eddington rates, whereas the high impact parameter simulations stay below 
Eddington accretion at all times. 



ing the accretion rates, from sub-Eddington accretion rates 
for large impact parameters to super-Eddington accretion 
for small impact parameters. The cloud velocity is the most 
important factor in determining the final Toomre unstable 



disk mass fraction, with a high mass for low velocity clouds 
and a low mass for high velocity clouds. 

Clearly, there is a lot of additional physics that must 
be included, which will be explored in a subsequent paper. 
Firstly, we did not investigate the thermodynamics (heat- 
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Figure 6. Logarithmic density in in the xy-plane with z=0 at different times for the test-case -^= 0. (d) is a close-up of the 

accretion flow visible in panel (c). The velocity is given in km/s. The whole cloud collapses into the black hole since the net angular 
momentum is zero. 



ing/cooling) of the disk as well as its fragmentation and star 
formation which is necessary to directly compare the sim- 
ulations with the observations. Secondly, a more realistic 
cloud model should be used since observations clearly indi- 
cate that GC clouds are highly turbulent ( Glist en fc Philippl 
I2OO4I ). Finally, black hole and stellar feedback will also play 
an important role in the formation process. 

Still our results are very encouraging, we have shown 
that redistribution of angular momentum must be taken into 
account when a cloud falls onto the GC black hole and our 
simulations show that this mechanism can indeed produce 
eccentric, Toomre unstable disks with properties similar to 
the ones observed in the GC. 

From all the simulations we have performed, simulation 
V03 shows the most promising features such as low eccen- 
tricity of around e — 0.24, a disk mass of roughly M — 10^ 
Mq, a size of r = 1.0 pc with sub-Eddington mass accretion 
at all times. 
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